Loss of purity by wave packet scattering at low energies 
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Abstract 

We study the quantum entanglement produced by a head-on collision between two gaussian wave 
packets in three-dimensional space. By deriving the two-particle wave function modified by s-wave 
scattering amplitudes, we obtain an approximate analytic expression of the purity of an individual 
particle. The loss of purity provides an indicator of the degree of entanglement. In the case the 
wave packets are narrow in momentum space, we show that the loss of purity is solely controlled 
by the ratio of the scattering cross section to the transverse area of the wave packets. 

PACS numbers: 03.67.Mn, 34.50.-s 



1 



In this Brief Report we describe the quantum entanglement generated by wave packet 
scattering in three-dimensional free space. Unlike one-dimensional problems studied pre- 
viously by one of us Q], scattering in 3D involves wave functions with much richer state 
structures for entanglement. Recently, we have demonstrated some interesting features for 
low energy eigen-functions in trapped systems Q]. For unbounded systems, Tai and Kurizki 
have analyzed the increase of entropy in terms of the scattering matrix|3|. Their approach is 
based on a particular form of two-particle wave functions in which the corresponding Schmidt 
decomposition can be expressed in pure plane wave bases |3j . For general two-particle wave 
functions, however, particles may not be paired in plane wave modes. Therefore a complete 
analysis of scattering effects on entanglement production remains open for investigations. 
Here we address the problem in the low-energy regime. Assuming the interaction potential 
is isotropic and short-ranged, we can employ the s-wave approximation to obtain the scat- 
tering wave functions. Our task is to determine the loss of purity of an individual particle, 
which serves as a measure of entanglement in our system with pure two-particle states. 

The system under investigation consists of two interacting particles of equal mass m in 
free space. The Hamiltonian in terms of center of mass and relative coordinates is given by: 
H = H cm + H rel with 

p2 

H cm = 2^ (1) 

Hrei = ^ + V ( r )- (2) 

Here M = 2m is the total mass and /i = m/2 is the reduced mass. For convenience, we will 
use the units with h = fi = 1. We assume that the interaction potential V (r) is isotropic 
and has a short range b such that V (r) w for r > b. Initially, the two particles are in the 
form of (disentangled) gaussian wave-packets, each having a width <To in momentum space. 
Their initial positions and average momenta are ±ro and =Fko respectively. The direction of 
k is chosen such that the packets make a head-on collision at later time (Fig. 1). 

The initial two-particle wave function in momentum space is given by a product state: 
$ (ki, k 2 , 0) = 0i (ki) 02 (k 2 ), where 



(ka) = r (k l5 k ; V2/a ) e-^-Xo)^ (3) 
(k 2 ) = T (k 2 , -k ; V2/a ) e ^+ k o) ? (4) 



with T (a, b; c) being a gaussian function parameterized by the inverse width c and the peak 




FIG. 1: An illustration of the system before and after a head-on collision of wave packets. Under 
the s-wave approximation, the scattered part of the single particle density (i.e., either particle 1 or 
particle 2) is a spherical shell shown in the right figure. The arrows indicate that the two particles 
go into opposite directions. 



at b, 

/c 2 \ 3/4 
r (a, b; c) = I — J exp 



(5) 



The function (5) allows us to express the wave packets in a compact form. After the 
scattering, the wave function in the long time limit takes the form: 

$ (ki, k 2 , t) = (N)- 1/2 [0! (kO 2 (k 2 ) e^+W + e<p scat (k 1; k 2 , t)} . (6) 

Here the first term corresponds to a non-scattering part that propagates freely, and the 
second term corresponds to the scattering part. The constants A" and e are normalization 
factors such that $ and 4> sca t are both normalized to unity. In this paper we treat |e| <C 1 

ClS db small number. 

To analyze the quantum entanglement, it is customary to study the entanglement en- 
tropy obtained from the Schmidt decomposition of (6). Quite generally, the Schmidt modes 
are not simply the momentum eigenfunctions, and the decomposition has to be performed 
numerically. We note that this is in contrast to the special case considered in Ref. [3], in 
which the Schmidt modes are momentum eigenfunctions. To gain insight of the problem 
analytically, we employ the purity function V as an alternative measure of entanglement. 
Such a function is denned by V = Tr(pf), where pi = Tr 2 (pi2) is the reduced density of 



the particle 1, and pi 2 corresponds to the two-particle density matrix associated with the 
state (6). For pure two-particle states considered in this paper, the smaller the value of V, 
the higher the entanglement. A disentangled (product) state corresponds to V = 1. We 
remark that V shares similar features as entropy, but it has the key advantage that it is 
more accessible to theoretical analysis 0,0,0]. In atomic physics, V (or its inverse V^ 1 ) 
las also been employed to indicate the two-body correlations in various dynamical processes 



Specifically, V takes an integral form in our system: 

V = J J I /^(ki'^Wka^^k!^ (7) 
From Eq. (4) and Eq. (6), we obtain the expression of P up to the second order of e: 



V*sl-2\e\*[l+I 1 -I 2 -I 3 ], 



(8) 



where the integrals ii, / 2 , h are defined by, 



h 



J J (jf scat (k l5 k 2 , t) 0! (kj) 2 (k 2 ) d'k^k, 
h = j J €cat (k 1; k 2 ,t)0 2 (k 2 )d 3 k 2 



d 3 k 



<fki. 



(9) 
(10) 

(11) 



These integrals describe the interference between a non-scattered wave and a scattered wave. 
It is interesting to note that there are no first order terms in e in Eq. (8), as these terms 
cancel each other once the e dependence in the normalization constant N is taken into 
account. We also remark that as long as e is a small parameter, Eq. (8) is valid for general 
two-particle states that are initially separable, not just for gaussian wave packets. 

To calculate e and <p sca ti let us rewrite Eq. (6) in terms of center of mass and relative co- 
ordinates: $ (ki, k 2 , t) = (N)~ 1/2 (p cm (K, t) 4> re i (k, t), where K = k x + k 2 , k = (ki - k 2 ) /2, 



(K,t)=r(K,O;l/0b) e 



Since we are interested in low-energy scattering processes, 



we may keep only the s- waves of the scattered part in cf) re i (k, £), i.e. 



J rel 



(M) 



iff 



(k,t)+e V ^ (k,t) 



(12) 



where 0^ e f (k,t) 



'rel 



(k,0)e 



is the freely propagating non-scattered part. In this 



way, we have </w(ki, k 2 , t) 



(K,t)^)(k,t). 



Under the assumption that the two particles are well separated (initially and finally) 
the scattering part er)^ (k, t) is given by, 

ev [s) (k, t) = -i-j / / <P rel (k', 0) / (AO — e^'"^ W^r 
(2n) j J r 



(13) 



where, 



/o(*) 



,i20(fc) _ ]_ 

2ik 



(14) 



is the s-wave scattering amplitude, and 9 (k) is the s-wave scattering phase shift. After some 
calculations, we obtain, 



ev (s) (M) 



(Tn 



4k - i2a^r 



2ie(k) 



r ( k, k k; — 



e^ fc - fc °) r ° _ ilk 2 t 
k e 2 ■ 



(15) 



The constant e is determined from the norm of the right side of Eq. (15). We may Taylor 
expand 9(k) at ko to the second order. The normalization condition is a Gaussian integral 
that can be calculated explicitly. This gives 

aid' (k ) 2 



kh 2 



1-Reie 



2i9(k ) 



cxp 



2-ial6"{h) 



(16) 



2fc 



in the 



V 2 - iVgfl" (ifeo 

where 9'(ko) and 9"(ko) are first and second derivatives of 9 (ho), and 7 2 = 1 
denominator corresponds to the spreading factor of the spatial width of the packets (since 
To / ko is the time of collision) . Therefore the spreading of the wave packets would decrease 
the norm of the scattered wave function \e\ 2 as expected. We also note that the value of the 
bracket [1 — Re {. . .}] is bounded between and 2, and therefore \e\ 2 is smaller than 2a 2 jk\. 

With the results of e and (j)^ s \ we find that I\ ~ ~ ~ ^ are all °f the 

order of ctq/Zcq. Because of the prefactor \e\ 2 in Eq. (8), these integral's contribution to P 
is about ctq//cq, which will be neglected. Hence the purity of final state is approximately 



1-2 lei 



(17) 



where \e\ 2 is given by Eq. (16). 

Further simplification of this result can be made in the limit al9" (ko) <C 1, and 
&o9' (ko) <C 1, i.e., the wave-packets are very narrow in the momentum space. In this 
limit, we have |e| ~ \^2\f (k )\ a c , where l/cr c = 7/er is the spatial width of the wave 
packets at the collision time. Alternatively, we may employ the scattering cross section 
So (ko) = 4vr|/ {ko) | 2 , so that 

ct 2 S (k ) 



l-V^4ai\fo(k ) 



7T 



(18) 



Therefore the purity of the two-particle wave function after scattering can now be explicitly 
expressed in terms of the s-wave scattering cross section as well as the widths of wave packets. 
However, we remark that such a simple relation is valid if <Jq is sufficiently small. The result 
can become more complicated when <jq9' (ko) or a\Q" (k ) in Eq. (16) are not negligible. 

Equations (16-18) are the main results of this paper. We see that the degree of entan- 
glement (quantified by 1 — V) is determined by a simple dimensionless parameter a^S(k ). 
Since l/a c is the spatial width of an individual wave packet at the collision time, a^S(k ) is 
just the ratio of scattering cross section to the characteristic cross area of the wave-packet 
in position space. Hence, a stronger entanglement can be generated for systems with a 
larger value of the ratio. For example, this can be achieved by exploiting resonance scatter- 
ing in which S(ko) can be enhanced near the resonance energies defined by the interaction 
potential |ll 0]. 

We point out that the degree of entanglement is typically small. This is due to the fact 
that the two-particle wave function is dominated by an un-scattered part, which is a product 
state. However, if mainly the scattered part is observed (for example, by detecting directions 
different from the incident one), then the relevant wave functions can have a much higher 
degree of entanglement. For the s-wave function given in Eq. (15), if o"o is small such that 
the phase shift can be treated as a constant 9(ko), then the normalized (relative coordinate) 
scattered wave function is a spherical shell of radius ko and thickness <7o in momentum space. 
We find that the corresponding purity function V has a leading term proportional to o\jk^ 
when ao/ko <C 1 is a small parameter. Therefore the narrower the width of the wave packet, 
the stronger the entanglement in the scattered part of the wave function. 

To conclude, we present a simple and general formula that approximates the loss of purity 
due to a head-on collision between two gaussian wave packets in three dimensional space. As 
long as the scattering is dominated by s-waves, our results provide a quantitative measure 
of quantum entanglement generated. In particular, our approach allows us to identify the 
key parameter a^S(ko), that explicitly connects the scattering cross section and the width 
of wave packets to the degree of quantum entanglement. 
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